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ABSTRACT 

We model early star forming regions and their chemical enrichment by Population III 
(Pop III) supernovae with nucleosynthetic yields featuring high [C/Fe] ratios and 
pair-instability supernova (PISN) signatures. We aim to test how well these chemical 
abundance signatures are preserved in the gas prior to forming the first long-lived 
low-mass stars (or second-generation stars). Our results show that second-generation 
stars can retain the nucleosynthetic signature of their Pop III progenitors, even in 
the presence of nucleosynthetically normal Pop III core-collapse supernovae. We find 
that carbon-enhanced metal-poor stars are likely second-generation stars that form in 
minihaloes. Furthermore, it is likely that the majority of Pop III supernovae produce 
high [C/Fe] yields. In contrast, metals ejected by a PISN are not concentrated in the 
first star forming haloes, which may explain the absence of observed PISN signatures in 
metal-poor stars. We also find that unique Pop III abundance signatures in the gas are 
quickly wiped out by the emergence of Pop II supernovae. We caution that the observed 
fractions of stars with Pop III signatures cannot be directly interpreted as the fraction 
of Pop III stars producing that signature. Such interpretations require modelling the 
metal enrichment process prior to the second-generation stars’ formation, including 
results from simulations of metal mixing. The full potential of stellar archaeology can 
likely be reached in ultra-faint dwarf galaxies, where the simple formation history may 
allow for straightforward identification of second-generation stars. 

Key words: dark ages, reionization, first stars - stars: Population II - stars: Popu¬ 
lation III - stars: chemically peculiar - galaxies: high-redshift - galaxies: formation 


1 INTRODUCTION 

The formation of the hrst, so-called Population III (Pop III), 
stars marks an important yet incompletely understood pe¬ 
riod in cosmic history. As the hrst non-linear baryonic struc¬ 
tures in the Universe, Pop III stars initiate the reionization 
and chemical enrichment of the early intergalactic medium, 
thus setting the stage for the assembly of more mature galax¬ 
ies (e.g. Bromm et al. 2009; Loeb & Furlanetto 2013, and 
references therein). 

The character of Pop III star formation has been theo¬ 
retically investigated in great detail, and simulations of this 
process have begun to converge on consistent conclusions 
(Bromm 2013; Glover 2013; Greif 2015). Unfortunately, ob¬ 
servations that can directly constrain these ideas are in short 
supply. The James Webb Space Telescope (JWST) may be 
able to observe supernovae (SNe) from the most massive 
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Pop III stars (Hummel et al. 2012; Pan et al. 2012; Whalen 
et al. 2013; Chen et al. 2014; de Souza et al. 2014). The lack 
of observed low-mass Pop III stars in the Milky Way can 
potentially constrain the low end of the Pop III initial mass 
function (IMF), a constraint that will gain in strength with 
increasing survey sizes (Hartwig et al. 2015). The ionizing 
effect of Pop III stars may be seen in epoch of reioniza¬ 
tion measurements (Greif & Bromm 2006), although this 
prospect has significantly weakened in light of the recent 
Planck determination of a lower optical depth to Thomson 
scattering (Planck Collaboration et al. 2015). Arguably, in 
the absence of such direct probes, the best available con¬ 
straints on Pop III star formation may come from chemical 
abundances in the most metal-poor stars, which preserve the 
chemical composition of their birth clouds in their stellar at¬ 
mospheres (e.g. Karlsson et al. 2013; Frebel & Norris 2015). 
The ideal objects for probing Pop III star formation are 
the first metal-enriched low-mass stars (or second-generation 
stars), which contain only metals from Pop III stars. 


© 2015 The Authors 



2 Alex P. Ji, A. Frebel, V. Bromm 


It is not known with any certainty what abundance 
signatures would clearly indicate enrichment primarily by 
Pop III stars, but intriguing empirical hints have been found. 
Some unique abundance patterns have been identified in 
metal-poor stars, including the high fraction of carbon- 
enhanced metal-poor (CEMP) stars at extremely low iron 
abundance (e.g. Norris et al. 2013; Placco et al. 2014), and 
a strong odd-even effect characteristic of pair-instability 
supernovae (PISNe) (Heger & Woosley 2002; Aoki et al. 
2014). However, average metal-poor halo stellar abundances 
([Fe/H] < —2.5) suggest that the majority of Pop III stars 
died as core-collapse supernovae (CCSNe) (Joggerst et al. 
2010). As the ejecta from multiple SNe mix together, any 
unique Pop III abundance patterns which could be preserved 
in metal-poor stars will be diluted away to an average core 
collapse SN yield. 

Recent simulations have suggested that second- 
generation stars should generically form by combining the 
metal yields from more than one Pop III SN explosion, even 
in the simplest sites of second generation star formation 
(Ritter et al. 2015; Jeon et al. 2015). Such an immediate en¬ 
richment by multiple SNe could rapidly dilute away unique 
Pop III chemical signatures in the majority of the subse¬ 
quent, Pop II, star formation events. Modeling such dilu¬ 
tion from multiple SNe is a classical problem in the theory 
of chemical evolution. Several approaches have been em¬ 
ployed, including models that assume instantaneous mix¬ 
ing (e.g. Tinsley 1980; Kirby et al. 2011), extensions that 
allow inhomogenous mixing (e.g. Oey 2000; Argast et al. 
2004; Karlsson 2005), and full hydrodynamic simulations 
(e.g. Kobayashi & Nakasato 2011; Webster et al. 2014). 
These models were all devised to study chemical evolution 
in a static self-enriching system over multiple generations, 
such as the Milky Way disc or a large dwarf galaxy, with 
the exception of Webster et al. 2014 who specifically inves¬ 
tigate chemical evolution in ultra-faint dwarf galaxies. How¬ 
ever, second-generation stars form in a unique initial stage of 
chemical enrichment, taking place in the earliest, low-mass 
building blocks of galaxy formation. The existing chemical 
evolution models focus less on the initial chemical enrich¬ 
ment stage and more on the subsequent chemical evolution. 

In this paper, we investigate how mixing the ejecta from 
multiple SNe may remove tell-tale nucleosynthetic signa¬ 
tures from Pop III stars. We consider second-generation star 
formation in the simplest possible environments motivated 
by cosmological simulations of the early Universe. Consid¬ 
ering these environments, we apply idealized mixing models 
to predict when the Pop III signature is preserved by metal- 
poor stars and can be probed by stellar archaeology. 

The structure of our paper is as follows. In Section 2, we 
discuss what is known about where and how stars form in the 
early Universe. In Section 3, we describe mixing models to 
produce chemical distributions both for second-generation 
stars and subsequent generations of stars. The results of 
our models are presented in Section 4, and we conclude in 
Section 5. Whenever needed, we use a cosmology with pa¬ 
rameters = 0.3, Ha = 0.7, h = 0.7, and Hj, = 0.05. 
The baryon fraction then is fb = 1/6. When gas is ionized 
and neutral, we use /^ionized = 0.59 and /^neutral = 1.23, re¬ 
spectively, where j-inip is the mean molecular weight of the 
gas. 



t (Myr) 



t (Myr) 

Figure 1. Schematic early star formation rates. Top: minihalo 
star formation rates. Pop III forms IOOMq in a 10® yr burst at 
virialization (2 = 25). Pop II forms IOOOMq in a 10 Myr burst 
after a recovery time of 20 Myr. Bottom: atomic cooling halo star 
formation rates. Pop III forms IOOOMq uniformly over 150 Myr 
following the merger tree in Greif et al. (2008). A Pop II burst at 
virialization (2 = 10) forms IOOOMq of stars over 10 Myr. 


2 STAR FORMATION IN THE EARLY 

UNIVERSE 

To model the first chemical enrichment events, we must 
consider the nature of star formation in the early Uni¬ 
verse. Here, we take an ab initio approach to understanding 
early star formation based on theoretical arguments and re¬ 
view the insights obtained from simulations. The picture 
described here motivates the mixing models presented in 
Section 3. 

Structure formation in ACDM is a hierarchical process, 
with small dark matter haloes forming early and merging 
into larger haloes. Two environments stand out in the early 
Universe where first star formation can take place. These en¬ 
vironments correspond to the two most significant cooling 
mechanisms in metal-poor gas. The first stars in the Uni¬ 
verse form in minihaloes, where star formation is triggered 
by molecular hydrogen cooling (e.g. Tegmark et al. 1997). 
The first galaxies in the Universe likely form in atomic cool¬ 
ing haloes, when hydrogen line cooling becomes significant 
(e.g. Oh & Haiman 2002). Cosmological simulations of these 
two environments have outlined a general picture for the for- 
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Figure 2. Top: distribution of number of Pop III SNe per mini¬ 
halo, assuming IOOMq of Pop III stars form. This was computed 
as described in the text from 10® minihaloes, assuming stars with 
M is IOM 0 turn into SNe. Bottom: Pop II star formation rates 
for a self-enriching atomic cooling halo. Model A is a 1 Gyr flat 
star formation rate. Model B is a 100 Myr burst. Model C is a 
500 Myr burst. 

mation of Pop III stars and the first few generations of Pop II 
stars, which we describe in Sections 2.1 and 2.2. The star 
formation rates described are shown schematically in Fig¬ 
ure 1. We discuss nucleosynthetic signatures from the hrst 
generations of stars in Section 2.3. 


2.1 Population III star formation 

2.1.1 Pop III initial mass function 

Pop III stars form in minihaloes of mass ~ IO^Mq around 
2 ; ~ 25, where star formation is triggered by molecular 
hydrogen cooling (Couchman & Rees 1986; Haiman et al. 
1996; Tegmark et al. 1997). The lack of metal cooling and 
the resulting higher temperatures suggest that Pop III stars 
have high characteristic masses, but our understanding of 
the characteristic mass has changed considerably over the 
last decade (e.g. Bromm 2013). Early simulations found that 
these stars might have characteristic masses > lOOM© with 
one star forming per minihalo (Abel et al. 2002; Bromm 
et al. 2002; Yoshida et al. 2006; O’Shea & Norman 2007). 
These stars would typically explode as PISNe with explo¬ 


sion energies > 10®^ erg (e.g. Heger & Woosley 2002). The 
strong feedback from these stars likely prevents further star 
formation within their host minihaloes (Greif et al. 2007; 
Whalen et al. 2008). 

However, more recent studies with higher spatial reso¬ 
lution have found that fragmentation may occur in the disc 
surrounding a Pop III protostar, forming small multiples 
of stars with lower characteristic masses > IOMq (Clark 
et al. 2008, 2011; Stacy et al. 2010; Greif et al. 2011, 2012). 
The amount of fragmentation varies between different mini¬ 
haloes and depends on the detailed thermal history of the 
gas (Greif et al. 2012). This general picture has been corrob¬ 
orated by observations of metal-poor stars, which seem to 
show clear a-element enhancements that are consistent with 
core-collapse supernovae (e.g. Beers & Christlieb 2005). 

Although a shape and range of the IMF cannot yet be 
determined robustly, Greif et al. (2011) find that the proto- 
stellar mass function across five minihaloes is approximately 
logarithmically flat ranging from ~ 0.1 — IOMq based on 
1000 years of fragmentation and accretion, although accre¬ 
tion should continue for 10® or 10® years in total. Susa et al. 
(2014) evolve 59 minihaloes with lower resolution for ~ 10® 
years after sink particle formation, also finding a fairly loga¬ 
rithmically fiat IMF ranging from ~ 1 —IOOMq. Such simula¬ 
tions suggest that the Pop III IMF is close to logarithmically 
flat from ~ 1 —IOOMq (with a characteristic mass ~ 2 OM 0 ), 
though a flat IMF is also reasonable (e.g. Cooke & Madau 
2014)^. With such an IMF, it is likely that multiple CCSNe 
form within each minihalo. Thus, we use a log-fiat IMF with 
lower and upper mass bounds of 1 and IOOMq. 

The mass range suggested here does not allow for higher 
mass PISNe. If PISNe occur, the PISN progenitor is likely 
the only Pop III star to form in a given minihalo. There is 
currently only tentative evidence for metal-poor stars with 
PISN abundance patterns (Aoki et al. 2014), because PISN 
signatures may only be found in higher metallicity stars and 
thus be missed by current metal-poor star surveys which 
preferentially select the most metal-poor stars (Karlsson 
et al. 2008). In addition, superluminous supernovae which 
might be associated with PISNe have been observed at both 
high and low redshift (Gal-Yam et al. 2009; Cooke et al. 
2012). If a PISN explodes in a minihalo, the violent feed¬ 
back associated with it will likely prevent any subsequent 
star formation there. 

2.1.2 Pop III star formation rate 

The Pop III star formation in an individual minihalo should 
occur in one burst. Collapse begins around the virialization 
time of the dark matter halo, and stars form in about a 
free fall time. The mass accretion rates into the central star 
forming regions are around 10“® —10 “^Mq yr“® for ~ 10® yr 
(e.g. Greif et al. 2011; Hirano et al. 2014; Susa et al. 2014). 
Not all of this mass ends up in protostars due to radiation 
feedback (Stacy et al. 2012), and each minihalo ends up 
forming around lOOM© of stars (e.g. Susa et al. 2014). 

Although our Pop III IMF allows low mass stars to 
form, for the investigation of metal enrichment we are only 

^ A logarithmically flat IMF has dN/dlog M oc constant, while a 
flat IMF has dN/dM oc constant. 
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interested in high mass stars which explode as SNe. The 
exact number of SNe in a minihalo is not well determined; 
simulations which resolve significant fragmentation do not 
run long enough to determine the actual stellar masses (e.g. 
Greif et al. 2011), while simulations running for longer times 
may not fully resolve the fragmentatien (e.g. Hirano et al. 
2014; Susa et al. 2014). As a simple way to model the number 
of Pop III SNe exploding within a minihalo, we sequentially 
draw stellar masses from the log-flat IMF ranging over 1 — 
IOOMq until we reach a total of M* = 50, 100, or 2 OOM 0 of 
stars. If adding a star would cause us to increase the mass 
beyond M*, we redraw from the IMF. Only stars with mass 
larger than lOM© are assumed to undergo a supernova, so 
we stop when we have drawn > M, — lOM© of stars. The 
distribution of the number of SNe found this way from 10® 
minihaloes is shown in Figure 2. 

This procedure gives initial stellar masses in addition to 
the total number of SNe in a minihalo, and in principle the 
stellar masses can be used as inputs to SNe nucleosynthesis 
calculations to determine the metal yields on a star-by-star 
basis. However, given the uncertainties in these calculations 
(see Section 2.3), we aim to separate the impact of SN yields 
and energies from the problem of metal mixing. Thus, we 
do not use the stellar masses derived from this procedure 
for anything other than estimating the number of SNe in a 
minihalo. 

Pop III stars will not form in atomic cooling haloes, 
unless in rare cases, because the haloes will typically be pol¬ 
luted with metals prior to their virialization (e.g. Johnson 
et al. 2008). One can instead estimate the number of star 
forming minihaloes that eventually merge into an atomic 
cooling halo. Greif et al. (2008) followed one atomic cooling 
halo and found that the ejecta from ~ 10 star forming mini¬ 
haloes mixed together during the formation of the atomic 
cooling halo, with the Pop III star formation occurring from 
t — 150 — 300 Myr. If each minihalo produces ~ IOOMq 
of Pop III stars, the Pop III star formation rate is about 
IOOOM 0 in 150 Myr. For an extreme upper bound, a maxi¬ 
mum of one hundred IO^Mq minihaloes can make up a single 
10 ®Mq atomic cooling halo. 

2.2 Population II star formation 

In this section we describe the formation of Pop II stars. The 
transition from Pop III to Pop II is often demarcated by a 
critical metallicity. The critical metallicity could be set by 
atomic line cooling (e.g. Bromm & Loeb 2003; Frebel et al. 
2007) or dust thermal cooling (e.g. Schneider et al. 2012; 
Ji et al. 2014). For simplicity, we do not apply any criti¬ 
cal metallicity, and instead assume all metal-enriched gas is 
able to form low mass Pop II stars that could correspond to 
metal-poor stars. In practice, our modelled metallicities are 
nearly always above Z ~ 10 , which places all gas above 

the critical metallicity for low-mass star formation according 
to the dust cooling criterion. 


2.2.1 Pop II star formation environment 

Both the minihalo and atomic cooling halo environments 
may host early Pop II star formation. Which environment is 
relevant in a given region of space depends on the degree of 


Table 1. Second-generation star forming environments 



Minihalo 

Atomic Cooling Halo 

Halo Mass 

10 ® Mq 

1O®M0 

^vir 

25 

10 

Rvir 

120 pc 

1.3 kpc 

'i^vir 

6 km 

18 km 

T ■ 

vir 

1900 K 

17400 K 

Dt^ [kpc^Myr“^] 

2.4 X 10-® 

8.1 X 10-4 


“ effective diffusion coefficient for turbulent mixing 


Pop III feedback and the corresponding recovery time (Jeon 
et al. 2014). 

If Pop III feedback is strong, then the minihalo loses 
most of its gas (e.g. Whalen et al. 2008) and the next time 
gas can collapse is in atomic cooling haloes with virial tem¬ 
perature ~ 10"^ K or mass ~ IO^Mq (Oh & Haiman 2002; 
Greif et al. 2008; Wise et al. 2012). These haloes are per¬ 
haps the best candidates for sites of first galaxies (Bromm 
& Yoshida 2011). If instead Pop III feedback is weak, then it 
is possible for gas to recollapse in the original minihalo and 
form Pop II stars there (e.g. Ritter et al. 2012, 2015; Cooke 
& Madau 2014). These stars form out of a mix between the 
supernova ejecta that have fallen back and primordial gas in¬ 
flowing to the centre of the minihalo through filaments. The 
minihalo may gain some mass during the recollapse (Jeon 
et al. 2014). 

It is likely that differing accretion/merger histories and 
the distribution of Pop III stellar masses result in differ¬ 
ent recovery times, allowing second-generation stars to form 
from both of these environments. The atomic cooling halo 
and minihalo will be the two main contributors to the pop¬ 
ulation of second-generation stars in the ACDM paradigm. 
However, there are two other situations that may arise. It is 
possible that the blast from a PISN in one minihalo is able 
to pollute a nearby minihalo and trigger second-generation 
star formation (e.g. Greif et al. 2010; Wise et al. 2012; Smith 
et al. 2015), although the degree of metal pollution may not 
be very high (Gen & Riquelme 2008). Also, star formation 
may be triggered in the shocked shell of a Pop HI super¬ 
nova explosion (e.g. Mackey et al. 2003; Chiaki et al. 2013). 
These two cases are possible variants on the picture of metal 
recollapse in a minihalo which depend on the specific hydro¬ 
dynamics of metal mixing and supernova shocks, and we do 
not expect the majority of second-generation stars to form 
through such mechanisms. However, such cases may be able 
to sample the metal output of individual SNe. 

As canonical values, we consider minihaloes of mass 
~ 1O®M0 virializing at 2 : ~ 25, and in case of moder¬ 
ate stellar feedback we assume second-generation (Pop H) 
star formation takes place in a similar environment. For the 
atomic cooling halo case, we assume a typical halo mass of 
~ 1O®M0, virializing at 2 ~ 10 (Greif et al. 2008). Note 
that in both these environments, there will be a distribu¬ 
tion of halo masses and virialization redshifts (Wise et al. 
2012; Hirano et al. 2014; Susa et al. 2014). The halo mass 
is especially important in the minihalo case as a higher halo 
mass increases the likelihood that second-generation stars 
can form within the minihalo (Cooke & Madau 2014; Jeon 
et al. 2014). Canonical physical parameters of these haloes 
are given in Table 1. 
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Regardless of whether a second-generation star forms 
in a minihalo or an atomic cooling halo, it is very likely 
that metals from multiple Pop III SN will contribute to 
the second-generation star forming gas, for the following 
reasons. If the second-generation star forms in an atomic 
cooling halo, its metal content is determined by mixing to¬ 
gether metals from many progenitor minihaloes. If instead 
the second-generation star forms in a minihalo, the Pop III 
stars which formed in that minihalo were of lower mass, 
resulting in less vigorous feedback. Thus, it is likely that a 
small multiple of Pop III stars contributed metals to second- 
generation stars. 

2.2.2 Pop II star formation rates 

When Pop II star formation occurs in a minihalo, the 
SN ejecta are mixed with pristine infalling gas, accreting 
through filaments, in a turbulent core, collapsing after a 
recovery time of ~ 20 Myr (Jeon et al. 2014; Ritter et al. 
2015). The rate of Pop II star formation in this case has not 
been simulated, but in the study by Ritter et al. (2012) the 
cool dense gas core at the centre has a mass ~ SOOOMq. If 
the star formation efficiency is ~ 30 per cent, then this will 
ultimately produce ~ lOOOM© of stars in ~ 10 Myr. After 
10 Myr, the more massive Pop II stars will explode as su¬ 
pernovae (~ 5SNe). If a sufficient number of SNe explode 
and the minihalo is not too massive, this will finish the gas 
expulsion initiated by the Pop III stars and halt star forma¬ 
tion in the minihalo until it merges into another system. On 
the other hand, if the minihalo has gained enough mass to 
contain the SN explosion, then the minihalo can continue to 
self-enrich and form stars out of gas contaminated by Pop II 
SNe (Jeon et al. 2014). 

Pop II star formation in an atomic cooling halo be¬ 
gins approximately when the halo virializes {z ~ 10 or 
t ~ 400 Myr) and atomic line cooling triggers star forma¬ 
tion (e.g. Greif et al. 2008; Safranek-Shrader et al. 2014b,a). 
The first burst of star formation will occur near the centre of 
the atomic cooling halo, after which star formation will be 
temporarily halted. Stars forming after this initial ~ 10 Myr 
burst of Pop II star formation will no longer be truly second- 
generation stars, having been polluted by the first generation 
of Pop II stars. The amount of stars formed in this burst is 
not certain. Greif et al. (2008) find that 10®M© of cold dense 
gas has collapsed with Z ~ 10~^Zq. Safranek-Shrader et al. 
(2014a) follow collapse to a gas density of 10^ cm“^ and find 
that stellar clusters can form with ~ lOOOM© of stars after 
4 Myr. Safranek-Shrader et al. (2015) resolve the same initial 
conditions down to a density of ~ 10^^ cm“®, resolving indi¬ 
vidual protostellar cores as sink particles and including the 
effect of their radiative feedback. In 20,000 yr, ~ SOM© of 
stars formed, although many of the protostars will continue 
to accrete gas beyond the simulated time. Thus it appears 
that approximately lOOOM© of second-generation stars can 
form in this initial 10 Myr burst. The SNe from this burst 
likely contaminate the dense gas with metals from Pop II 
stars (Wada & Venkatesan 2003), preventing further second- 
generation stars from forming afterwards. 

After the initial burst, Pop II star formation in the 
atomic cooling halo will continue until it is halted. This 
could happen externally through reionization or merging 
into another halo, or it could happen internally through SN 


Table 2. CEMP signature supernova yields 



CCSN“ 

Faint SN® 

Wind-|-SN= 

Pop II SN'* 

Mass (M©) 

20 

20 

20 

20 

®SN,51® 

1 

0.74 

1 

1 

C (M©) 

0.211 

0.20 

1.034 

0.128 

Fe (M©) 

0.072 

1.09 xlO-® 

0.072 

0.073 


“ Nomoto et al. 2006 
Iwamoto et al. 2005 

CCSN with added C from Hirschi 2007 (see text) 

Nomoto et al. 2013 

® i?SiM 51 is the explosion energy in units of 10®^ erg 

feedback. We estimate these processes with three models. 
Model A is a constant star formation rate for 1 Gyr followed 
by an abrupt cutoff. Model B is a short ~ 100 Myr burst. 
Model C is a longer ~ 500 Myr burst. The burst in models 
B and G use the equation 

V'W = i’o (1) 

where ipo is some normalization and r is the characteris¬ 
tic time of the burst (100 Myr and 500 Myr for models B 
and C, respectively). These star formation rates are plot¬ 
ted in Figure 2. The normalizations are all set to produce 
30,000 M©kpc"®. 

2.2.3 Pop II initial mass function 

Modern star formation theory has recognized the impor¬ 
tance of supersonic turbulence in regulating star formation 
(Mac Low & Klessen 2004; McKee & ©striker 2007). The 
turbulent velocity helps prevent global gravitational col¬ 
lapse, while the complex network of shocks formed by super¬ 
sonic flows create local density fluctuations that can trigger 
star formation. The turbulent density spectrum may play a 
fundamental role in determining the IMF (Hopkins 2012). 

In atomic cooling haloes, supersonic turbulence is cre¬ 
ated during virialization of the halo, and during subse¬ 
quent accretion from the cosmic web (e.g. Greif et al. 
2008; Safranek-Shrader et al. 2014a,b). This suggests that 
the Pop II IMF in atomic cooling haloes may be similar 
to that in present-day star formation. We thus take the 
IMF in atomic cooling haloes to be the Salpeter IMF from 
0.1 — 100 M©. To date, there have not been significant stud¬ 
ies of turbulence in minihaloes which have experienced rec¬ 
ollapsing gas. We assume the IMF in this situation is also 
Salpeter. 

2.3 Nucleosynthetic signatures 

We now discuss the Pop III chemical abundance signatures 
that may be found in second-generation stars, as well as 
other early metal sources. Unfortunately, calculating super¬ 
novae yields from first principles is a difficult task. Core col¬ 
lapse supernovae in particular have proven to be difficult to 
understand from first principles due to a complex explosion 
mechanism. After the star’s core becomes gravitationally un¬ 
stable and undergoes core bounce, the outgoing shock stalls 
and is presumably revived by some combination of hydro- 
dynamic instabilities and neutrino heating (e.g. Janka et al. 
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Table 3. PISN signature supernova yields 



CCSN“ 

PISN*' 

Pop 11 SN" 

Prog. Mass (M 0 ) 

20 

195 

20 

Esn,51 

1 

40 

1 

C (M 0 ) 

0.211 

4.13 

0.128 

Na (M 0 ) 

0.0029 

0.00028 

0.00181 

Mg (M 0 ) 

0.150 

4.39 

0.247 

Ca (Mq) 

0.00623 

0.993 

0.00921 

Fe (M 0 ) 

0.072 

3.08 

0.073 

Co (Mq) 

1.5 xlO-'^ 

5.59 xlO-® 

6.21 X 10“® 

Ni (Mq) 

0.00175 

0.00825 

7.11 X 10-* 


“ Nomoto et al. 2006 
Heger & Woosley 2002 
Nomoto et al. 2013 


2012). 3D simulations of SN explosions that treat multifre¬ 
quency neutrino transport are a recent development, and the 
results of these calculations have not yet converged between 
different groups (e.g. Hanke et al. 2013; Mezzacappa et al. 
2014). Instead of a full hydrodynamic simulation, CCSN 
metal yield calculations often parametrize the explosion in 
different ways (e.g. Kobayashi et al. 2006; Heger & Woosley 
2010; Limongi & Chieffi 2012). These methods agree quali¬ 
tatively, but not in detail (Nomoto et al. 2013). 

Rather than looking to reproduce detailed elemental 
abundances, we instead focus on the most important broad 
signatures of Pop III nucleosynthesis. Many Pop III CCSN 
calculations tend to produce abundances consistent with the 
stellar halo chemical record (e.g. Tumlinson 2006; Joggerst 
et al. 2010), suggesting that Pop III CCSNe produce a stan¬ 
dard a-enhanced yield not qualitatively different from their 
higher metallicity counterparts. As a representative yield, we 
take the 2 OM 0 , zero-metallicity model with a SN explosion 
energy of Esn = 10®^ erg from Nomoto et al. (2006). 

However, there are two abundance signatures that may 
distinguish themselves from this typical yield pattern, the 
carbon enhanced metal-poor (CEMP) signature and the 
pair instability supernova (PISN) signature. We now discuss 
these signatures in more detail. 

The exact yields we use in our future models are given 
in Tables 2 and 3. We have quoted the yields with a high 
degree of precision, but the yields will likely span a range of 
values. Our adopted values are only meant to be represen¬ 
tative realizations. 


2.3.1 CEMP signature 

If a star has [C/Fe] > 0.7, we say it displays the CEMP sig¬ 
nature following Aoki et al. (2007)^. It is not known whether 
the CEMP signature is created by Pop HI stars, or if there 
is a non-primordial explanation related to surface pollution 
or a critical metallicity (see Norris et al. 2013, and Nomoto 
et al. 2013, for a comprehensive account of possibilities, and 
Sluder et al. 2015 for a mechanism not previously consid¬ 
ered). We proceed under the assumption the CEMP signa¬ 
ture does indeed trace nucleosynthetic yields of Pop III stars. 
There are two classes of mechanisms to create the signature. 

2 [X/Y] = logj^giNx/Ny) - logio(A''x/Ny )0 for element X,Y 


It can be produced either by removing iron from or adding 
carbon to the standard CCSN yields. 

A likely mechanism for removing iron is to allow some 
sort of metal fallback during the SN explosion. This can be 
accomplished through faint supernovae that undergo mixing 
and fallback (Umeda & Nomoto 2002, 2003; Iwamoto et al. 
2005) or a jet-like explosion (Tominaga et al. 2007). As a 
representative fallback yield, we take a faint CCSN with 
Esn = 0.74 X 10®^ erg from Iwamoto et al. (2005) which re¬ 
produces the detailed abundances of HE 1327—2326 (Frebel 
et al. 2005). 

A likely mechanism for adding carbon to the yield is by 
winds from rotating massive Pop III stars (Meynet et al. 
2006, 2010; Hirschi 2007). To approximate this, we add 
O.823M0 of carbon to our standard Pop III CCSN. This 
represents the carbon winds from a very fast rotating star 
of 2 OM 0 (Hirschi 2007), and results in [C/Fe] = 0.89 for the 
wind model. 

2.3.2 PISN signature 

The PISN signature is still a theoretical signature, although 
it has been tentatively observed in a metal-poor star (Aoki 
et al. 2014). In contrast with CCSNe, the explosion mecha¬ 
nism for PISNe is better understood and does not need to be 
put in artificially (e.g. Heger & Woosley 2002). The typical 
PISN turns a large fraction of its stellar mass into metals 
and exhibits suppressed odd-even abundance ratios. We will 
consider a star to have the PISN signature if [Co/Ni] < —0.5, 
where Co and Ni are representative odd and even elements. 
We take the yields from Heger & Woosley (2002) for a 
I 95 M 0 Esn ~ 40 X 10®^ erg explosion as a characteristic 
PISN yield. 

2.3.3 Other metal sources 

It will be of interest to look at how Pop II SNe dilute away 
the CEMP and PISN abundance signatures. For Pop H SN 
yields, we use the 2 OM 0 , Z = 0.001 yields from Nomoto 
et al. (2013) which are similar to the Pop HI CCSNe yields. 

There are two other main sources of metals in standard 
chemical evolution models: AGB stars and Type la SNe. 
Low-to-intermediate mass stars undergo an AGB phase, dur¬ 
ing which they expel significant amounts of carbon in their 
winds. A typical time until an intermediate-mass star un¬ 
dergoes the AGB phase is ~ 100 Myr, much longer than the 
recovery time if second-generation stars form in minihaloes, 
but shorter than the recovery time for an atomic cooling 
halo. However, the total metal mass from a single AGB 
star’s winds is not very large, and its importance in standard 
chemical evolution is largely due to the steep slope of the 
IMF. The top-heavy Pop HI IMF suggests that AGB stars 
will not be an important contributor to elements found in 
second-generation stars. Type la SNe produce large amounts 
of iron-peak elements. In principle, there could be Pop HI 
binary white dwarfs that produce Type las. These also are 
unlikely to contribute to second-generation stars, since such 
systems may be rare, and the delay time between their for¬ 
mation and detonation is typically ~ 1 Gyr (although it 
could be as low as 100Myr, see Maoz et al. 2012). For these 
reasons, we will not consider contributions from these two 
metal sources. 
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Figure 3. Illustrating the enrichment model for an atomic cool¬ 
ing halo (Section 3.2). The Mixed Lagrangian Volume (MLV) has 
a Lagrangian radius Rj^ = 6.4 kpc at 2 = 10. Not shown is the 
size of the ACH Lagrangian Volume (ALV), which has a radius 
^L.ACH = 11.4 kpc at 2 = 10. Here, three minihaloes have expe¬ 
rienced SN explosions. In this example, q\ = 0, 92 = 1, and 53 
is between 0 and 1. The distances di are uniformly drawn from 
a spherical volume V (not shown) whose radius is larger than 
Rl + max{Ri}. 


3 METAL ENRICHMENT MODELS 

As seen in Section 2, most second-generation stars form 
ont of gas that has been enriched by multiple SNe, even 
in the simplest possible environments of minihaloes (MH) 
and atomic cooling haloes (ACH). Thns, there are two fun¬ 
damental questions we must answer to understand when a 
unique Pop III chemical abundance signature can be pre¬ 
served: (1) How much does the initial mixing of multiple 
Pop HI SN ejecta dilute this signature? And, (2) how long 
does it take for Pop H SNe to erase the Pop III chemical 
signature? 

We answer these questions with the following frame¬ 
work. Consider two types of SNe, a normal type n and a 
special type s. Let ps be the probability that a given SN is 
special. Then the yields from SN i are given by: 

f Vs probability 

Vi = \ (2) 

[yn probability 1 — Pa, 

where yi is the vector of metal yields for SN i in units of solar 
masses {not the mass fraction). In this paper, y„ will always 
denote the metal yields from standard Pop HI CCSNe (i.e. 
the first column of Tables 2 and 3). ys will be the metal 
yields for faint SNe, wind+SNe, or PISNe. 

To determine the fraction of stars preserving the special 
signature during the initial enrichment period, we run these 
yields through an enrichment model that determines the 
metallicity of second-generation stars. We consider a mass 
Afmix within which the gas is chemically homogeneous due 
to turbulent mixing. If the jith second-generation star is en¬ 
riched by Nj Pop HI SNe, then the metallicity Zj is given 
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by 


2 .= 




Nj 

i=0 


( 3 ) 


where qi G [0,1] is the fraction of metals from the ith SN 
that is mixed into Mmix. Depending on the environment, 
there will be different ways to determine Mmix, Nj, and Qi. 
We describe simple models for these quantities for second- 
generation star formation in minihaloes (Section 3.1) and 
atomic cooling haloes (Section 3.2). 

To examine question (1), we simplify the problem as 
follows. For a given metallicity Zj, we classify the star as 
special or normal depending on its metal ratios. When in¬ 
vestigating the CEMP signature, we consider the star to be 
special if [C/Fe] > 0.7. When investigating the PISN signa¬ 
ture, we consider the star to be special if [Co/Ni] < —0.5. For 
each value of Ps, we employ a Monte Carlo procedure to run 
our enrichment model Wotai times to produce {Zj}^^’^'^'. 
Each of these Monte Carlo realizations represents a chemi¬ 
cally homogeneous star cluster, of which Ns star clusters are 
classified as special. Then we define the preservation fraction 
fa to be the fraction of these star clusters that retain the 
special signature: 


fs{Ps) = 


Ns 

Ntotal ’ 


( 4 ) 


where we have made explicit the dependence of fa on the in¬ 
put probability for special yields, pa- The function fa{Pa) in 
different environments thus describes how well a special sig¬ 
nature is preserved, depending on the likelihood for special 
SNe to occur. In this paper we use Wotai = 10® throughout. 

To answer question (2), in Section 3.3, we describe when 
the contribution from Pop H SNe is able to erase the unique 
Pop HI signature. 


3.1 Second-generation star formation in 
minihaloes 

All Pop HI SNe in a minihalo explode near the centre, driv¬ 
ing a cumulative blastwave. If fallback occurs and triggers 
second-generation star formation, some fraction of the met¬ 
als from each SN will enrich this second generation (Ritter 
et al. 2012; Jeon et al. 2014). When modeling the chemi¬ 
cal abundances of metal poor stars, the most common as¬ 
sumption is that metals from different SNe stay completely 
unmixed, and abundance ratios reflect yields from a single 
SN (e.g. Joggerst et al. 2010; Nomoto et al. 2013). This ide¬ 
alized case implies fa{Pa) = Pa- However, as we argued in 
Section 2, this is likely not the generic case, even in mini¬ 
haloes. In general, ~ 1 — 10 SNe will enrich a minihalo. We 
select Nj from the distributions in Figure 2. 

To mix multiple SNe, the simplest assumption is that 
all metals from all SNe completely mix into some gas mass 
(e.g. Tinsley 1980; Cooke & Madau 2014). In the language 
of Equation 3, this fully-mixed model implies qi = 1 for all 
SNe. 

However, a more interesting and likely case is that not 
all metals from each SN will fall back. Ritter et al. (2015) in¬ 
vestigate the scenario of recollapse in a minihalo after seven 
SNe have exploded in the centre. About half the SN ejecta 
escape from the halo, suggesting a similar fraction of the 
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E (erg) 



r (kpc) 


Figure 4. Top: SN bubble final radius expanding out of a mini¬ 
halo as a function of initial energy deposited and ionizing feed¬ 
back. The ionizing feedback is parametrized in terms of ngat, 
which is the flat density in the inner part of the champagne solu¬ 
tion. The dotted line indicates the virial radius of the minihalo. 
Bottom: density profiles for different riflat- The temperature is 
10"^ K throughout, which affects the gas pressure. 


metal mass escaping. The recollapsing gas in these simnla- 
tions is biased towards metals from the first SNe to explode 
in the minihalo, since this ejecta runs into a denser surround¬ 
ing medium and thus is able to cool faster. The difference in 
fallback is a factor of 2 — 3 between the first and last star to 
explode. The amount of metal fallback depends on the de¬ 
tails of star formation, stellar feedback, and hydrodynamics 
in each minihalo. This suggests that a stochastic model is 
required to describe the amount of fallback. For simplicity, 
we model qi in the partial metal fallback case as a random 
variable uniformly distributed between 0 and 1. To approx¬ 
imate the effect of different SN energies, we also consider a 
biased case where all the metals from faint SNe are retained 
in the minihalo {qi = 1 if the ith SN is special), but the nor¬ 
mal SNe only retain a uniformly distributed fraction of their 
metals {qi ~ Unif[0,1] if the ith SN is normal). A suite 
of more detailed calculations of metal mixing and fallback 
such as those suggested in Ritter et al. (2015) may find that 
other fallback distributions are more realistic. 


3.2 Second-generation star formation in atomic 
cooling haloes 


Metal enrichment for second-generation star formation in 
an ACH is a complex process that generally requires a full 
hydrodynamic simulation (e.g. Greif et al. 2010; Wise et al. 
2012; Jeon et al. 2015). However, given the significant pa¬ 
rameter uncertainties in Pop III star formation, we here 
describe a much simpler model that captures the physical 
nature of this process. The model is illustrated in Figure 3. 

Consider an ACH of virial mass Mach ~ IO^Mq that 
virializes at Zvir = 10. If we trace all particles within the 
ACH to higher redshift, there is an associated Lagrangian 
volume that encompasses all the mass within the ACH at 
Zvir- We call this volume the ACH Lagrangian Volume, or 
ALV for short. This Lagrangian volume is conceptually the 
same Lagrangian volume that is sampled at high-resolution 
in a zoom-in simulation (e.g. Safranek-Shrader et al. 2014b). 
If the Lagrangian volume is spherical, there is some La¬ 
grangian radius Rl.ach such that 


Rl.ach 


/ 3Mach \ 

V 47rpm ) 


(5) 


where pm is the total matter density evaluated at some red- 
shift. The physical size of this Lagrangian radius at «vir is 
found by using Pm(zvir), giving Rl.ach = 11.4kpc. 

Now consider the volume V surrounding and including 
the ALV. We assume that star-forming minihaloes form uni¬ 
formly distributed within V, the SNe in these minihaloes ex¬ 
plode independently, and the spherical SN remnant bubbles 
do not interact with each other hydrodynamically. Given the 
density of star-forming minihaloes n, the number of mini¬ 
haloes in the volume V is given by a Poisson distribution 
with mean nV. We draw Nj from this Poisson distribution. 

Turbulence in the ACH will uniformly mix metals in 
much of the ACH gas (Greif et al. 2010), but not all the 
gas within the ACH is well-mixed (i.e. Mmix SS /tMAcn). 
We show in Section 3.2.2 that Mmix ~ 3 x lO^M©. The 
Lagrangian volume corresponding to the mixed gas is cor¬ 
respondingly smaller. We call this volume the Mixed La¬ 
grangian Volume, or MLV. The radius of the MLV is denoted 
Rl, a smaller radius than Rl,ach- 

For a single minihalo, the fraction of its metal content 
that is mixed into the MLV can be estimated by calculat¬ 
ing the volume overlap between its SN bubble and the La¬ 
grangian volume of the MLV. In other words, if the ith SN 
bubble has radius Ri, the MLV has radius Rl, and the cen¬ 
tre of the minihalo is located at a distance di from the centre 
of the MLV, then 


Ho 

47r/3Rf 


( 6 ) 


where the overlap volume Ho can be found by simple geom¬ 
etry and is given by 


° ^ ^ - rzdi - 

r rjS 

-Ui 


di> R 

sum. 

RdifI < di < Rs um ? 

di ^ Rdiff 

(7) 


where Rsum = Rl + Ri and Rdis = \Rl — Ri\- For each 
minihalo, we draw di from a uniform 3D spatial distribu¬ 
tion, where the maximum distance is given by taking V as a 
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sphere. A uniform distribution is appropriate as the distance 
is considered prior to the collapse of the ACH (see Figure 3). 
Although the SN bubble radius expands with time, we do 
not consider the time dependence and instead let the ra¬ 
dius instantaneously expand to its final (fixed) radius, only 
depending on the SN explosion energy. This is further ex¬ 
plained in Section 3.2.1. 

Other than pa, the only parameters needed are n and 
V. We choose n so that Amh, the mean number of mini¬ 
haloes whose centre is in the ALV, is 10, 30, and 100. This 
corresponds to n « 0.0016, 0.0048, and 0.016 minihaloes per 
kpc® at z = 10. We choose V to be a sphere around the MLV 
whose radius is 3.5 kpc larger than the MLV (8 kpc when 
we increase the SN bubble radius beyond that). Since the 
number of SNe is Poisson, the exact choice of V does not 
matter as long as it is large enough to enclose Rl + Ri ■ 

In practice, the mass of the SN bubble is likely concen¬ 
trated in a thin shell, and mixing into the spherical volume 
depends on the development of Rayleigh-Taylor instabilities 
(Madau et al. 2001). If no such instabilities develop and all 
the metals are concentrated in the shell, the overlap fraction 
is instead determined by the overlapping solid angle, given 
by: 


fo 

J 2Ridi-df-Iif+Rl 

1 4R^ 

u 


di> R 

sum 

Rdifl < di < R sum 
di ^ Rdiif 


( 8 ) 


We have checked that using the thin-shell formula does not 
make a significant difference in our results. 

Since multiple SNe explode in each minihalo, for each 
SN bubble we combine the yields from ~ 1 — 10 SNe with 
the M* = IOOMq number distribution in Figure 2. 


3.2.1 Supernova bubble radius 

To calculate the size evolution of a Pop III SN bubble, we 
use the thin-shell approximation in an NFW halo (Madau 
et al. 2001). We consider haloes of mass IO^Mq virializ- 
ing at z = 25. The NFW concentration parameter is deter¬ 
mined using the model from Correa et al. (2015) to be 2.6, 
although note that minihaloes are strongly triaxial and a 
spherical NFW density may not properly characterize the 
gravitational potential (Sasaki et al. 2014). Different con¬ 
centrations do not greatly affect the bubble radii. 

To properly model the gas density profile, we must ac¬ 
count for ionizing radiation from the Pop III progenitor stars 
(Bromm et al. 2003). To approximate the typical gas den¬ 
sity profile resulting from ionizing radiation, we use the self¬ 
similar champagne flow solution of Shu et al. (2002) as im¬ 
plemented by Wang et al. (2012). We parametrize the self¬ 
similar scale using ngat, the density in the inner region. A 
larger nuat corresponds to a smaller stellar lifetime (Alvarez 
et al. 2006). Outside of the shock, we use an isothermal 
profile (p oc r”^), evaluated at the halo virial temperature. 
When the isothermal halo density reaches the ambient IGM 
value, we prevent it from decreasing further. We use a gas 
temperature of lO'^ K throughout, corresponding to ionized 
gas since the HII region outpaces the SN blast wave. We 
use the cooling function of Sutherland & Dopita (1993) for 
collisional ionisation equilibrium with [Fe/H] = —3.0. The 
bubble is initialized with a Sedov-Taylor solution at 10'* yr. 



r (Myr) 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

M(10® Mq) 


Figure 5. Top: mixing mass within an atomic cooling halo. Solid 
lines indicate the mixing due to turbulent diffusion, which we use 
to determine the MLV. Light gray line indicates our choice for 
the MLV of 3 X 10® Mq. Dashed lines indicate how the mixing 
mass changes if including a SN bubble of 200 pc, which is rele¬ 
vant for calculating the Karlsson (2005) model. Bottom: mixing 
mass distributions for the three Pop II star formation rates in 
the bottom panel of Figure 2. Solid, dashed, and dotted lines cor¬ 
respond to models A, B, and C respectively, k is the number of 
Pop II SNe enriching a given low mass star. All distributions are 
roughly between 1 — 2 X IO^Mq. 


Unlike the original Madau et al. (2001) model, we include 
only one SN burst, so we do not include the time-dependent 
mechanical luminosity term. Instead, SNe contribute their 
energy to the initial energy of the SN bubble. We take the 
stopping radius to be the point where the shock velocity 
is equal to the sound speed in the photoionized medium 
(T = 10"* K). Although this model is simple, it reproduces 
the more detailed calculations of Greif et al. (2007), where 
the SN bubble stalls at 2.5 kpc for a 10®^ erg explosion in a 
5 X 1O®M0 minihalo. 

The results are shown in Figure 4. The top panel shows 
the SN bubble radius as a function of total energy in the bub¬ 
ble, while the bottom panel shows the gas density profiles 
that the bubble pushes against. Not shown is the underly¬ 
ing NFW profile, which slows the bubble gravitationally. If 
a bubble is energetic enough to break past the champagne 
flow shock, it is able to travel much further into the IGM 
before stopping. In all cases with E ^ 10®^'® erg, the bubbles 


MNRAS 000, 1-18 (2015) 










10 Alex P. Ji, A. Frebel, V. Bromm 


reach their maximum radius in less than 100 Myr. Since the 
time delay between Pop III star formation and virialization 
of the ACH is ~ 100 Myr (see Figure 1), this justifies using 
an instantaneous bubble size in our ACH enrichment model. 
Based on Figure 4, a good fiducial radius for a 10®^ erg 
CCSN explosion is ~ 0.2 — 0.4 kpc, while a 10®^ erg PISN 
explosion reaches ~ 3 kpc. We use 0.2 kpc and 3 kpc as the 
inital radius for CCSNe and PISNe, respectively. However, 
if the SN bubble expands far enough from the minihalo to 
be swept up by the ambient Hubble expansion, the radius 
can grow by a factor of up to 2.5. The maximum extent of 
the SN bubbles will then be Ikpc and 7.5 kpc for CCSNe 
and PISNe, respectively. We show models using these larger 
radii as well. 


3.2.2 Mixing mass 

The size of the Mixed Lagrangian Volume in an atomic cool¬ 
ing halo is primarily determined by the extent of turbulent 
diffusion (Greif et al. 2010). To model this, we use the mix¬ 
ing volume formula from Karlsson et al. (2008) 

Hmix(r) = ^ (6Ar)®/' (9) 

where r is the time over which turbulent diffusion occurs, 
and Dt is the turbulent diffusion coefficient. The factor of 
6 is a combination of a factor of 2 from the ID solution to 
the diffusion equation and a factor of 3 from three spatial 
dimensions (Karlsson et al. 2013). The turbulent diffusion 
coefficient is given as 

Dt ~ (Uturb)4urb/3 (10) 

where Uturb is a typical turbulent driving velocity, and fturb 
is a typical turbulent driving scale. This is an order-of- 
magnitude estimate based on dimensional grounds, simi¬ 
lar in nature to that used by mixing length theory to de¬ 
scribe convective transport in stars where the mixing length 
is parametrized in terms of a characteristic scale height (e.g. 
Kippenhahn et al. 2012). Since turbulence is driven by grav¬ 
ity in our haloes (Wise & Abel 2007; Greif et al. 2008; 
Safranek-Shrader et al. 2012; Ritter et al. 2015), we esti¬ 
mate that the characteristic Uturb is given by the halo virial 
velocity, and the characteristic £turb is given by a fraction of 
the virial radius ~ i?vir/10. Values for the diffusion coeffi¬ 
cient are given in Table 1. The mixing mass as a function of 
T and the ambient gas density are shown in Figure 5. 

A typical gas density in an atomic cooling halo is 
n ~ 10“^cm~®, and a typical amount of time for the gas 
to be turbulently mixed is the halo free fall time, tff ~ 
0.1/H{z) « 70 Myr for 2: = 10. This results in a typical 
mixing mass of Mmix = 3 x lO^M© of gas, or Rl = 6.4 kpc. 

Note that there is in principle a lower bound to the size 
of the MLV, since a single star cluster is chemically homoge¬ 
nous (e.g. Bland-Hawthorn et al. 2010; Feng & Krumholz 
2014). In the Milky Way this is emprically found to be 
~ 1O®M0, but it could be lower in the earliest galaxies 
(Safranek-Shrader et al. 2014b). 


second generation stars from low mass stars that have been 
contaminated by Pop H SNe. This can happen in massive 
minihaloes as well as in atomic cooling haloes, where po¬ 
tential wells are sufficiently deep to retain their gas if SN 
feedback is not too strong. At this point, the halo is massive 
enough to constitute a self-enriching system, a characteristic 
of a first galaxy (Bromm & Yoshida 2011). To answer ques¬ 
tion (2), how quickly Pop H SN can erase unique Pop HI 
signatures, we add Pop H SNe one at a time until the signa¬ 
ture is wiped out. Here, as before, it is important to know 
how much gas the SN ejecta mixes into. 

For a massive minihalo, in our spherically symmetric 
model a single 10®^ SN fills the entire virial radius. Thus 
we can instantly mix the Pop II SNe into the entire mini¬ 
halo volume. For simplicity, we assume the total gas mass 
available is the same as in the original minihalo, although 
in principle gas will be continually accreted. 

In an atomic cooling halo, the mixing mass for an in¬ 
dividual SN is often assumed to be the maximum radius 
of a supernova remnant (e.g. Ryan et al. 1996). However, 
turbulent mixing can expand the mixing mass. To estimate 
the size of this effect, we use the Karlsson (2005) model to 
calculate mixing mass distributions using the three Pop H 
star formation models in Figure 2. We employ the mixing 
volume formula from Karlsson et al. (2008), which is similar 
to Equation 9 but now accounts for the inital SN bubble size 
Rsn- 

Knix(r) = ^ {6DtT + Rl^f^ . (11) 

We use Rsn = 200 pc from Figure 4, and n = 0.12 cm“® and 
Dt = 7X 10~* kpc^ Myr“^ from Karlsson (2005). The results 
for a small number An of Pop H SNe are plotted in Figure 5. 
Increasing Nu allows for slightly higher mixing masses, but 
the typical mixing masses are 1 — 2 x 1O®M0, only slightly 
larger than the mixing mass right after the SN bubble stops 
expanding. Lowering the star formation rate by a factor of 
10 only slightly increases the maximum mixing mass here 
to ~ 3 X 1O®M0. Thus, it is reasonable to assume that the 
mixing volume for Pop H SNe is about the mass of its SN 
remnant, which we take to be 3 x 10®M0 to maximize the 
dilution. It is only at larger Nu that the effect of turbulent 
mixing will be important. 

We can also come to the conclusion that turbulent mix¬ 
ing is relatively unimportant by estimating the time-scale for 
turbulent mixing to dominate the mixing volume in Equa¬ 
tion 11. This is approximately when 


^SN 

ODt 


27 Myr 


Rsn 
200 pc 


Dt 


10“^ kpc^Myr 


( 12 ) 

This time-scale is longer or comparable to the ~ 10 Myr 
recovery time in a self-enriching system, so turbulent mixing 
will not affect the mixing mass by more than a factor of a 
few. 


4 RESULTS 


3.3 Contamination by Pop II supernovae 

Since Pop H CCSNe and standard Pop HI CCSNe produce 
similar metal yields, it may be difficult to distinguish true 


We now present model calculations showing how well the 
CEMP (Section 4.1) and PISN (Section 4.2) chemical signa¬ 
tures are preserved in gas from which metal-poor stars would 
form in the early Universe. We inject a natal fraction ps of 


MNRAS 000, 1-18 (2015) 



Preserving primordial chemieal signatures 11 



Pfaint Pwind + CCSN 



Pfaiiit Pwind + CCSN 

Figure 6. CEMP signature preservation in star-forming haloes after the addition of metals by Pop III SNe. Top row: Minihalo models 
are shown. The solid, dashed, and dot-dashed lines indicate different amounts of SN ejecta fallback into the; minihalo. Different colours 
indicate different total mass of Pop III stars formed per minihalo, which translates into the SN number distributions in Figure 2. Overall, 
the observed /cbmp in second-generation stars reflects the input special fractions ps- Bottom row: atomic cooling halo models are shown. 
The solid and dashed lines indicate different SN bubble sizes. Different colors indicate different values for the mean number of minihaloes 
contributing metals to an ACH. In both cases, /cbmp rapidly changes from 0 to 1 over a narrow range of input ps- See text for more 
discussion. 


special Pop III yields along with 1 — Ps of normal Pop III 
CCSN. We process this through our mixing models (see Sec¬ 
tion 3), and calculate the fraction fa of stars retaining the 
special Pop III signature. We consider the preservation both 
during the initial enrichment by combining multiple Pop III 
SNe, as well as after adding Pop II SNe. In Section 4.3, we 
compare our models from Section 4.1 to the observed stellar 
halo CEMP star fraction. 


4.1 Preserving a CEMP signature 


As a reminder, we say a metal-poor star displays the CEMP 
signature If It has [C/Fe] > 0.7. The input special fraction ps 
will refer to two different models for producing the CEMP 
signature in Pop III stars. If pa is pfaint, then we combine 
faint and normal CCSNe to produce the CEMP signature. 
If Pa is Pwind-i-ccsN, then we add a carbon wind to the yields 
of some of the CCSNe to produce the CEMP signature. 
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Figure 7. CEMP signature preservation after dilution by A^u Pop II SNe. Panels are arranged as in Figure 6. Different line thickness 
and transparency denote different numbers of Pop II SNe. The MH lines show deviations from the model using M* = IOOMq of Pop III 
stars per minihalo and full metal fallback in the minihalo after adding Pop II SNe. The ACH models show deviations from the model 
using Rsn = 200 pc and IVmh = 30. In all models, only a few Pop II SNe are required to eliminate the CEMP signature, even if 100 per 
cent of Pop III stars produce the CEMP signature. Note that in several panels, the larger A^n lines are at /cemp ~ 0 regardless of ps 
and thus are barely visible at the bottom of the panels. 


4-1.1 Multiple Pop III SNe 

In Figure 6, we show the preservation fraction for the CEMP 
signature in the MH and ACH environments (top and bot¬ 
tom rows, respectively) assuming the signature is generated 
through either the faint-|-normal CCSN or the wind-|-CCSN 
mechanisms (left and right columns, respectively). As an 
example of interpreting the figure, consider the model of 
faint-bnormal CCSNe in a MH (top-left corner). Suppose 
Pfaint = 0.6, i.e. 60 per cent of all Pop III SNe are faint 
CCSNe and 40 per cent are normal CCSNe. Then the solid 
green line shows that if minihaloes form a combined IOOMq 


of Pop HI stars and all the metals fully mix (qi = 1), then 
~ 35 per cent of metal-poor stars forming in minihaloes will 
retain the CEMP signature. For reference, the thin black 
line in the plots indicates where /cemp = Ps, or where the 
CEMP signature has been perfectly preserved. We stress 
that all the second-generation stars in any given minihalo 
will have the same signature, since they form from a chemi¬ 
cally homogeneous star cluster (Ritter et al. 2015; Safranek- 
Shrader et al. 2015). 

In the top row of Figure 6 we show the CEMP signa¬ 
ture preservation fraction for the minihalo models. In gen¬ 
eral, the MH Wind+CCSN models are within ~ 10 per cent 
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of perfect preservation. However, the faint+normal CCSNe 
models can deviate strongly from perfect preservation. In 
particular, at low pfaint and larger numbers of Pop III SNe, 
the CEMP signature is strongly suppressed. This is simply 
because it is easier to reduce the [C/Fe] ratio by adding iron 
from the normal CCSNe to the iron-poor faint CCSNe. As 
expected, having biased fallback for faint SNe (dot-dashed 
lines) results in much better signature preservation. 

In the bottom row, we show the CEMP signature preser¬ 
vation fraction for the atomic cooling halo models. The solid 
and dashed lines show SN expansion radii of 200 pc and 
1 kpc, showing the range of preservation fractions based on 
whether the Hubble flow expands the SN bubbles. The col¬ 
ors indicate different numbers of minihaloes that contributed 
metals to the ACH. Recall that each minihalo in these mod¬ 
els can produce multiple CCSNe, where the number is drawn 
from the green lOOM© distribution in the top panel of Fig¬ 
ure 2. While the MH models reasonably preserve the Pop HI 
signature regardless of pa, the ACH models show a sharp 
transition from poorly preserved to over-preserved. This is 
due to the average [C/Fe] ratio in all haloes being deter¬ 
mined by Pfaint- As we increase the number of SNe contribut¬ 
ing metals to a gas reservoir (MH or ACH), the abundance 
scatter between different gas reservoirs decreases, resulting 
in a larger fraction of the gas reservoirs being close to the 
average value. Note that in the IVmh = 10 case, when pfaint 
and Pwind+CCSN = 1 there is still not perfect preservation. 
This occurs because of cases where all minihaloes are within 
the ALV but do not intersect the MLV. 

Overall, it is generally the case that the CEMP sig¬ 
nature can be preserved in the gas after combining yields 
from multiple Pop HI SNe in MHs. It is difficult but not im¬ 
possible to preserve the signature in ACHs. However, even 
if second-generation low-mass stars can be identified from 
observations of metal-poor stars, our models show that the 
observed CEMP star fraction will generally be slightly lower 
than the actual fraction of Pop III SNe producing the CEMP 
signature. A metal enrichment model (such as the one used 
to calculate Figure 6) is required to convert observed CEMP 
star fractions into the actual fraction of Pop HI SNe produc¬ 
ing the CEMP signature. 

4-1.2 Adding Pop II SNe 

In Figure 7, we show how the preservation signature is kept 
after adding Pop H SNe yields to the metal distributions 
from Figure 6. Each Pop II SN yield is diluted into 3 x 
1O®M0, which we have conservatively chosen to be on the 
high end of possible dilution masses. Interestingly, in all four 
cases, adding 1 — 3 Pop H SNe is enough to erase the CEMP 
signature in the gas as imprinted by Pop HI SNe. In the MH 
case, this is because the signature is already close to being 
erased by a few normal Pop HI CCSNe. In the ACH case, it 
is a result of the ~ 10 times discrepancy in the mixing mass 
between Pop III and Pop H SNe. 

A typical Pop II burst in these environments produces 
~ IOOOMq of stars. For a Salpeter IMF, there are ~ 5 mas¬ 
sive stars that explode as SNe in the first Pop H star for¬ 
mation burst. As a consequence, this strongly suggests that 
if CEMP stars trace Pop III nucleosynthesis, they are truly 
second-generation stars that form out of only Pop HI metals 
(although multiple Pop III stars will contribute metals). 



PPISN 



PpiSN 


Figure 8. PISN signature preservation in the ACH model. Top: 
Pop III only. The line colors indicate different values for the mean 
number of minihaloes. Solid, dashed, and dot-dashed lines indi¬ 
cate different SN bubble radii. The thin black line denotes perfect 
preservation. If a minihalo exeriences a PISN, that is the only SN 
in the minihalo. If it has CCSNe, then the number of CCSNe is 
drawn from the IOOMq distribution of Figure 2. Note that ppisN 
only extends from 0.0—0.2, as we do not expect most minihaloes 
to host PISNe. Bottom: after adding a single Pop II SN, the PISN 
signature is erased. The Nn = 1 line has /pisn ~ 0 for all shown 
PpiSN and is thus barely visible at the bottom of the plot. 
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Median [Fe/H] Median [Fe/H] 

-2.88 -3.01 -3.18 -3.48 -3.48 -2.88 -2.88 -2.88 -2.88 -2.88 -2.88 



0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 

Pfaint Pvr'md + CCSN 


Median [Fe/H] Median [Fe/H] 

-3.61 -3.71 -3.84 -4.04 -4.33 -3.61 -3.61 -3.61 -3.61 -3.61 -3.61 



Pfaint Fwind + CCSN 

Figure 9. CEMP signature preservation compared to expected CEMP star fractions based on model [Fe/H] and observed CEMP star 
fractions (Equation 13). Shaded regions show the 25 and 75 percentile of the model [Fe/H] distributions translated into observed CEMP 
star fractions, with light solid line showing the median value. The median [Fe/H] for the IOOMq and IVmh = 30 models are shown on the 
top axis of each panel. Panels are arranged as in Figure 6. For the MH panels, we only show the full fallback case. For the ACH panels, 
we only show Usn = 0.2 kpc. The values of Pfaint a-nd Pwind+CCSN where the model lines intersect the shaded bars are consistent with 
the observed CEMP fraction. 


4.2 Preserving a PISN signature 

We define a metal-poor star as retaining the PISN signature 
if it forms from gas with [Co/Ni] < —0.5. As previously dis¬ 
cussed, this is a representative odd-even element abundance 
ratio. In the top panel of Figure 8, we show the preserva¬ 
tion fraction for the PISN signature (/pisn) depending on 
the fraction of minihaloes hosting PISNe (ppisn). If a mini¬ 
halo does not host a PISN, it then has a random number 
of CCSNe drawn from the IOOMq distribution of Figure 2. 
Models are shown for different SN bubble radii and different 


numbers of minihaloes. We only use the ACH model and not 
the MH model, since the strong PISN explosion evacuates 
its minihalo of gas (Whalen et al. 2008). 

In general, the PISN signature is well preserved for the 
expected number of 10 minihaloes contributing to the ACH. 
As more minihaloes are added, the number of normal Pop HI 
CCSNe increases and is able to overwhelm the PISN signa¬ 
ture at these low ppisN fractions. Note that when the size 
of the PISN bubbles is increased, the PISN signature is ac¬ 
tually stronger than perfect preservation. This is because 
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for the same number density of minihaloes (as parametrized 
by Nmh), increasing the SN bubble radius allows more such 
bubbles to overlap with the MLV. The MLV is then able 
to sample more PISN bubbles than CCSN bubbles. Increas¬ 
ing the size of the CCSN bubbles somewhat counteracts the 
PISN signature for the same reasons. 

The PISN signature is very fragile. The bottom panel 
of Figure 8 shows that a single Pop II SN is able to wipe out 
the PISN signature. Again, this is due to the vastly different 
dilution factors between the Pop III and Pop II SNe. 

It is interesting to compare the results of this model 
to cosmological simulations, such as the one by Jeon et al. 
(2015). In contrast with many previous simulations, this sim¬ 
ulation allows for both PISNe and CCSNe to occur in the 
minihaloes. By 2; « 10.5, 2 PISNe have formed out of a total 
of ~ 200 Pop III stars, corresponding to ppisN ~ 0.01. These 
two PISNe contribute the majority of metals in the simula¬ 
tion, but nearly all of these metals are ejected into the IGM 
rather than found in virialized haloes. Their most massive 
galaxy is thus predominantly enriched by CCSNe. We also 
find that the Nmh = 30 case (similar to their number of 
minihaloes) has CCSNe dominating the metal enrichment 
in a typical ACH when ppisN = 0.01. 

4.3 Interpreting observations of metal-poor stars 

Ultimately, the goal of stellar archaeology is to interpret 
abundances of metal-poor stars to learn about properties of 
Pop III stars. We now present an example of how our results 
in Section 4.1 can be used to interpret observed CEMP star 
fractions in the Milky Way halo as constraints on Pop III 
properties. 

Let us suppose that the CEMP star fraction in the stel¬ 
lar halo is representative of second-generation stars. In prac¬ 
tice, the halo is an amalgam of stars from many sources, and 
even at fairly low metallicities some halo stars may form af¬ 
ter the first few generations. Regardless, with this assump¬ 
tion, we can compare the observed CEMP star fraction to 
that obtained by our metal enrichment models. Our models 
give estimates of the [Fe/H] distributions of Pop II stars. We 
take these [Fe/H] distributions and find the 25th and 75th 
percentile values to define a range of /cbmp. We consider 
models to be allowed if their predicted /cemp is consistent 
with the /cemp computed from the model [Fe/H] distribu¬ 
tions. For the observed CEMP star fraction, we take an ap¬ 
proximate fit to the CEMP fraction from figure 16 of Placco 
et al. (2014) for [C/Fe] > 0.7. 

( 0.2 [Fe/H] > -2.5 

/cemp = < —0.8 — 0.4 X [Fe/H] —4.5 < [Fe/H] < —2.5 
[1.0 [Fe/H] <-4.5 

(13) 

In Figure 9 we show the preservation fraction for some of 
the models in Figure 6. The shaded bars in the background 
indicate the observed CEMP star fraction evaluated at the 
modelled [Fe/H] distributions. Note that the shaded bars 
for the wind-bCCSN models are flat, since every SN has the 
same iron yield. 

When the model lines overlap the shaded bars, the cor¬ 
responding pfaint or Pwind+CCSN are consistent with the ob¬ 
servations. Taking the figure at face value, we see that the 
MH models produce the expected carbon fractions, although 


it generally requires a fairly high ps- The faint-l-normal CC¬ 
SNe models require pfaint ^ 0.6. If this picture is correct, 
then our results imply that > 60 per cent of Pop HI stars 
exploded as a faint SN. The high fraction of faint SNe can 
be interpreted as an indirect tracer of the IMF (e.g. Heger & 
Woosley 2002). A similar calculation for the biased fallback 
model suggests that Pfaint ^ 0.5. If instead the CEMP signa¬ 
ture is produced in MHs by carbon winds, the wind-|-CCSN 
models require 0.7 > Pwind-i-ccsN > 0.4. This could be in¬ 
terpreted as the fraction of Pop HI stars that are rotating 
rapidly enough to have significant winds. 

In contrast, the ACH models seem to be unable to pro¬ 
duce consistent CEMP star fractions without a huge num¬ 
ber of minihaloes contributing to the ACH. Recall that 
each minihalo here can produce multiple Pop HI SNe, and 
Nmh = 100 is a strong upper limit assuming there has been 
no smooth mass accretion. Thus, if the CEMP signature is 
created solely by Pop HI stars, a large fraction of Pop HI 
stars must produce the high [C/Fe] abundances, and most 
of these stars form in minihaloes. 

Given the simplicity of our enrichment models, we cau¬ 
tion that more sophisticated models could somewhat affect 
the exact numbers. In particular, while the overall metal- 
licity in our models is likely correct to within an order of 
magnitude, the actual value can be shifted by changing the 
dilution mass within some reasonable range. From Equa¬ 
tion 13, a 0.2 dex change in the model [Fe/H] results in a 
0.08 change in the observed /cemp- R is also a fairly strong 
assumption that the observed CEMP stars represent only 
second-generation stars, although at [Fe/H] < —3 it is prob¬ 
ably reasonably accurate. However, we can qualitatively see 
the effect of including subsequent generations of stars. Since 
the CEMP signature is rapidly erased by Pop H SNe, any 
contaminating stars will reduce the observed CEMP star 
fraction. Thus a true second-generation CEMP star fraction 
will be higher than the shaded bars plotted in Figure 9. This 
further increases the Ps required to be consistent with the 
MH models, as well as exacerbating the inability of the ACH 
models to produce the observed CEMP star fractions. Re¬ 
gardless, this approach of comparing observed special frac¬ 
tions fs to those found from models is an emerging frame¬ 
work to test where our understanding needs to most improve 
to constrain properties of Pop HI stars. 


5 SUMMARY AND CONCLUSIONS 

In this study we focused on whether unique Pop HI chem¬ 
ical signatures can be preserved in the gas from which the 
first low mass stars form. Using simple models that cap¬ 
ture the essence of initial Pop HI chemical enrichment, we 
calculate the fraction of second-generation stars preserving 
unique Pop HI chemical signatures depending on the input 
fraction of Pop HI stars that produce said signature. We 
now summarize our main findings and conclusions. 

The chemical signature of Pop III can be pre¬ 
served in gas from which the first low-mass stars 
formed. Metal-poor stars are believed to trace the nucle- 
osynthetic yields of their progenitor stars which in the case of 
the most metal-poor ones are individual Pop III SNe events. 
While this is the fundamental assumption that stellar ar¬ 
chaeology rests on, detailed investigations to what extent 
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SN yield signatures remain intact in the gas after the explo¬ 
sion and factoring in various environmental aspects such as 
mixing and the presence of multiple SNe are still lacking 

Even the very first enrichment step has irreducible com¬ 
plexity, with at least a handful of Pop III SNe contribut¬ 
ing metals to most second-generation star forming environ¬ 
ments. We find that the mixing and re-accretion of gas onto 
the halo following these SNe does not always wipe out the 
original nucleosynthetic signatures, and thus that it is pos¬ 
sible for second-generation stars to trace Pop III properties. 
This validates the endeavor of stellar archaeology. Neverthe¬ 
less, it becomes increasingly clear that interpreting observed 
abundance signatures in the most metal-poor stars requires 
using chemical enrichment models to factor in effects such 
as metal mixing and turbulence (see also Ritter et al. 2015). 
Otherwise the properties of Pop III stars may not be cor¬ 
rectly recovered. 

The frequency of metal-poor stars with a spe¬ 
cial abundance signature (e.g. high [C/Fe]) cannot 
be directly interpreted as the fraction of Pop III 
stars producing that signature. By comparing the re¬ 
sults of our MH and ACH models to the observed fraction of 
CEMP stars in the Galactic stellar halo, we find that CEMP 
stars likely formed from gas in a minihalo enriched by a SN 
that underwent metal fallback. The observed high CEMP 
star fraction is not necessarily produced by faint SNe alone, 
but if so then the fraction of Pop III stars producing the 
CEMP signature must be rather high, at least ~ 40 — 60 
per cent even though the observed CEMP star fractions are 
lower. This high fraction could be an indication of a biased 
scenario, such as preferential second-generation star forma¬ 
tion in minihaloes that experience low energy SNe (Cooke & 
Madau 2014). However, a better understanding of the CCSN 
explosion mechanism is likely necessary to confirm this, as 
even quite energetic SNe can produce a similar abundance 
pattern (Tominaga et al. 2007). Either way, it is nearly im¬ 
possible to reproduce the observed CEMP star fraction in 
atomic cooling haloes. 

Any Pop III chemical signatures are quickly 
erased after the emergence of Pop II SNe. This fun¬ 
damental result suggests that stars with unique and un¬ 
usual chemical abundance signatures are clean probes of 
Pop III SN yields. This effect has been understood before 
in the general context of chemical evolution (e.g. Audouze 
& Silk 1995), but the sheer rapidity with which unique 
abundance signatures are erased has not been appreciated. 
Future abundance studies based on high-resolution spec¬ 
troscopy should focus on any unusual halo stars as they 
will likely reflect Pop III nucleosynthesis yields. Moreover, 
stars with [Fe/H] < —4.0 of which many show unique abun¬ 
dances can all be regarded second-generation stars. This is 
in line with other suggestions from SN yield fitting of ob¬ 
served abundance signatures in the most metal-poor stars 
(Tominaga et al. 2014; Keller et al. 2014; Placco et al. 2015), 
although these authors also considered stars with halo-like 
abundances and not just unusual signatures. 

A PISN signature is quickly wiped out by the 
emergence of Pop II SNe. A strong odd-even effect is 
believed to be a tell-tale sign of a PISN signature. So far, 
only one star has been identified that displays some degree of 
an odd-even effect (Aoki et al. 2014). This may suggest that 
PISNe are exceedingly rare. However, this does not neces¬ 


sarily imply that Pop HI PISN explosions were equally rare, 
or even absent. Our results instead suggest that any PISN 
enrichment pattern would rapidly be hidden under subse¬ 
quent layers of Pop H core-collapse enrichment. In addition, 
a PISN signature might only occur in stars with high metal- 
licity since the Ca yield of a PISN is rather large (Karlsson 
et al. 2008). Current selection biases on the basis of a weak 
Ca H K line would thus preclude finding these stars. Interest¬ 
ingly, even though PISNe produce large amounts of metals, 
it appears to be very difficult for most of the PISN metals 
to be retained in a first galaxy environment that would then 
lead to the formation of metal-poor stars. This is due to met¬ 
als being dispersed over a much larger volume compared to 
those from CCSNe. As previously suggested, the best place 
to search for the PISN signature may be in IGM gas at high 
redshifts with absorption spectroscopy (Jeon et al. 2015). 

Despite the fact that our enrichment models are ideal¬ 
ized and only capture the essence of the problem, our results 
are encouraging and highly relevant to the interpretation of 
the abundance signatures of metal-poor stars and thus to the 
field of stellar archaeology. Hydrodynamic simulations that 
elucidate the nature of metal mixing and enrichment in the 
early universe are particularly important. The results from 
these simulations will allow for more quantitative investi¬ 
gations since our simple models have already demonstrated 
that one can recover valuable information about Pop HI 
stars from stellar archaeology. 

We have shown that it is possible to learn about the 
nature and properties of Pop HI stars by using chemical 
abundance measurements of metal-poor halo stars, such as 
in the case of the observed CEMP star fraction. However, the 
complex formation history of the stellar halo makes it dif¬ 
ficult to isolate true second-generation stars. The only way 
to best choose candidates is to select the lowest metallicity 
stars, but that may not yield a complete sample since not 
all second-generation stars may be at the lowest metallicities 
(see e.g. Frebel & Bromm 2012; Frebel et al. 2014). 

A more promising way to find and study second- 
generation stars may be in dwarf galaxies and to carry out 
dwarf galaxy archaeology (Frebel & Bromm 2012). As each 
dwarf galaxy is a single, confined object with one specific his¬ 
tory, it will be easier to model the galaxies and to separate 
second-generation stars from other stars within these galax¬ 
ies. Of the Milky Way’s dwarf galaxies, the ultra-faint dwarfs 
may be the best places to clearly identify second-generation 
stars. Their chemical abundances and star formation histo¬ 
ries suggest that ultra-faint dwarfs form stars only over a 
very short period of time (Frebel et al. 2014; Brown et al. 
2014), and due to their inefficient star formation, a large 
fraction of their stars (or even all of them) may be second- 
generation. If ultra-faint dwarf galaxies were the descendents 
of a single Pop III star-forming minihalo, each ultra-faint 
dwarf galaxy would only preserve one single combination of 
Pop HI yields. However, if these galaxies are ACHs that form 
from a merger of more than one second-generation star form¬ 
ing minihalo, then multiple second-generation gas reservoirs 
may be represented in the dwarf galaxy. 

A disadvantage of using stars in dwarf galaxies to con¬ 
strain properties of Pop HI stars is the relatively small num¬ 
ber of stars with available chemical abundance observations. 
High resolution spectroscopy is the gold standard, but the 
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faintness of most dwarf galaxy stars may require further 
development of medium resolution techniques (e.g. Vargas 
et al. 2013, 2014; Kirby et al. 2015) or await the era of ex¬ 
tremely large telescopes. 

All but the least-luminous ultra-faint dwarfs have ex¬ 
perienced some level of extended star formation (Vargas 
et al. 2013). Thus, for detailed studies of true Pop III sig¬ 
natures, it will be important to separate bona-fide second- 
generation stars from subsequent star formation. A promis¬ 
ing approach may be to combine different Pop III yields and 
enrichment models with chemical evolution models that ac¬ 
count for hierarchical galaxy formation (e.g. Salvadori et al. 
2015). These models will need to account for the effects that 
recent hydrodynamic simulations have highlighted (Ritter 
et al. 2015; Sluder et al. 2015; Bland-Hawthorn et al. 2015). 
With a larger sample of available dwarf galaxy abundances 
and more sophisticated modeling of the formation of these 
systems, it may become possible to robustly and comprehen¬ 
sively constrain the properties of Pop III stars with dwarf 
galaxy archaeology. 
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